Assisted finite-rate adiabatic passage across a quantum critical point: 
Exact solution for the quantum Ising model 
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The dynamics of a quantum phase transition is inextricably woven with the formation of excitations, as a 
result of the critical slowing down in the neighborhood of the critical point. We design a transitionless quantum 
driving through a quantum critical point that allows one to access the ground state of the broken-symmetry phase 
by a finite-rate quench of the control parameter. The method is illustrated in the one-dimensional quantum Ising 
model in a transverse field. Driving through the critical point is assisted by an auxiliary Hamiltonian, for which 
the interplay between the range of the interaction and the modes where excitations are suppressed is elucidated. 
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The complexity involved in describing a generic many- 
body quantum system prompted Feynman to suggest the use 
of a highly controllable quantum system as a simulator of an- 
other, generally complicated, quantum system of interest [I]. 
From this perspective, interesting quantum systems are those 
with a large amount of entanglement and hardly tractable in 
classical computers [2]. Quantum simulation has become an 
exciting field of research, which is being developed experi- 
mentally by exploring a variety of platforms including ultra- 
cold atoms, trapped ions, photonic quantum systems and su- 
perconducting circuits, among others. Simulation of many- 
body interacting systems is particularly advanced in imple- 
mentations with trapped ions [ 3] where the building blocks of 
a digital quantum simulator for both closed [4] and open [5] 
quantum systems have been demonstrated. Moreover, while 
early experimental efforts have been limited to somewhat low 
number of qubits, the simulation of few-hundreds of spins 
with variable-range spin-spin Ising-type interactions has re- 
cently been reported {«\. 

In a continuous quantum phase transition, divergence of 
length and time scales across a quantum critical point (QCP) 
leads inevitably to non-adiabatic dynamics. When a param- 
eter A of the Hamiltonian is changed across its critical value 
Ac, the energy gap between the ground and the first excited 
state vanishes, and adiabaticity breaks down. The Kibble- 
Zurek mechanism (KZM), originally developed for classical 
and continuous phase transitions [7, 8], predicts that the re- 
sulting density of excitations obeys a power-law scaling with 
the quench rate. The power-law exponent is expressed using 
the critical exponents at equilibrium and the dimensionality 
of the system [0, 10]. As a result, quantum quenches are use- 
ful to characterize universal features of a system, and to shed 
some light on its dynamics out of equilibrium. 

The inevitable formation of excitation is however undesir- 
able for a wide range of applications, such as the preparation 
of novel quantum phases in quantum simulation, and adiabatic 
quantum computation. Suppressing excitations is also of in- 
terest to variety of operations in the laboratory, like entangling 
strings of atoms [11]. This has motivated studies including the 
use of the energy gap arising from the finite size of the sys- 



tem [12], optimal non-linear passage across a QCP [13, 14], 
inhomogeneous quenches [15-17], and optimal quantum con- 
trol strategies [18]. All those approaches can be regarded as 
strategies to exploit or engineer a spectral gap. Notwithstand- 
ing, there is a need for new methods to ensure adiabaticity [2]. 
In this letter, we shall exploit recent advances in the simula- 
tion of coherent A;-body interactions [5, 1 Q] and transitionless 
quantum driving [20, 21] to explore an alternative to quan- 
tum adiabatic protocols, and assist a fully adiabatic finite-rate 
passage across a QCP. 

Shortcut to the adiabatic driving of a two-level system. - 
Demirplak and Rice [20], and Berry [21] have shown the pos- 
sibility of implementing a transitionless quantum driving in 
multilevel systems. Let us consider the Landau-Zener (LZ) 
transition, the simplest model supporting the KZM [22], de- 
scribed by a Hamiltonian: 
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= A(f)cT' + Ac7-'^ 



(1) 



where cj'^ ^^^ are the usual Pauli matrices. The instantaneous 
eigenbasis reads: 

1 1 (A)) = sin 1 1 (-oc)) - cos |2(-oc)) , 
|2(A)) =cos0|l(-oo))+sin0|2(-oc)), 

where the angle Q obeys the relations 

A . A 



cos20 = 



VA2+A2 



sin20 — 



VA2+A2 



and the energy gap is £2(0 - Ex (f ) = 2VA2 + A2. Following 
[20, 21], it is found that the Hamiltonian that drives the exact 
evolution of the system along the adiabatic solution associated 
with the instantaneous eigenbasis {|«(A))} of //q in Eq. (1) is 
given hy H = Hq + Hi with [23] 



Hi=iX'it)Y,[\df.n){n\-{n\d^n)\n){n\]. 
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For the LZ crossing, upon explicit calculation one finds 

1 A 



Hi 



^X'{t] 



2A2 + A(f)2 



(2) 



(3) 
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The adiabatic solution of Hq, in which the instantaneous 
eigenstates pick up exclusively a phase along the evolution- 
the sum of the dynamical and Beny phases - becomes the ex- 
act solution of the time-dependent Schrodinger equation as- 
sociated with H ~ Hq + H\ no matter how fast the transition 
is crossed, i.e. how large the rate X'{t) is. This approach has 
recently been verified in the laboratory with an effective two 
state model arising in a Bose-Einstein condensate in the pres- 
ence of an optical lattice [24]. 

Models.- We turn now our attention to the family of d- 
dimensional free-fermion Hamiltonians: 



■^ = Lv^k[«k(A(f))-ak]V4, 

k 



(4) 



where Ok = ((7^,C7j^,a^) denote the Pauli matrices acting on 
the k-mode and !//][ = i'^ki^'^ki) fermionic operators. 
The function flk (A ) = (fl-^(A),fl-j^(A),flj^(A)) is specific for the 
model and the sum goes over independent k-modes. Such a 
Hamiltonian represents a variety of systems with QCP, includ- 
ing in particular the Ising and XY models ind=l [25], as well 
as the Kitaev model ind = 2 [26] and d = 1 [27]. As such it 
has been the subject of a recent series of works on defect pro- 
duction induced by a quantum quench [9]. 

We shall use it to illustrate and investigate the pos- 
sibility of driving an adiabatic passage across a QCP. 
Let us consider the instantaneous eigenstates of with 
eigenenergies associated with the k-mode ek,± — ±|ak('^)| ~ 



±y'a^(A)2+4(A)2 + 4(A)2. We generalize Eq. (3) to find 
the modified Hamiltonian = Mq + Jfi, where 



0^ 



k 



2e2 



[(flk(A) X dxdxJyk)) ■ o-k] i/k 



induces the adiabatic crossing of the QCP by driving the dy- 
namics exactly along the instantaneous eigenmodes of 

Without further knowledge of the explicit form of the ma- 
trix elements of J^, its form in real space cannot be deter- 
mined. Next, we turn our attention to a specific model. 

The quantum Ising model in a transverse field.- Consider a 
chain of N spins described by the \d quantum Ising model in 
a transverse magnetic field g. 



(5) 



a paradigmatic model to study quantum phase transitions [25] 
of relevance to current experimental efforts in quantum simu- 
lation [28]. We assume periodic boundary conditions On+i = 
(Ti and, for simplicity, even A^. This model exhibits a quan- 
tum phase transition at gc — ±\ between a paramagnetic phase 
{\g\ > 1) and ferromagnetic phase {\g\ < 1). 

The Jordan- Wigner transformation, = 1 — 2clc„, ((7,f + 
iOn) = 2c„n/<n(l — '^c]ci), where c„ are fermionic annihi- 
lation operators, allows us to rewrite the Hamiltonian (5) as 
a free fermion model. Below, we will limit ourselves to the 
plus-one-parity subspace of the Hilbert space, that includes 



the ground state - note that commutes with the parity op- 
erator P = n^Li <^n- In fermionic representation: 



n=\ 



with anti-periodic boundary conditions c^r+i = —c\. Us- 
ing the Fourier transform c„ ~ e I^Q-e'^'W^ with mo- 
menta consistent with the boundary conditions k E ~ 
{±7z/N,±3n/N,...,±{N — l)7z/N), we can conveniently 
rewrite it as 

= 2 ^ [alig - cosk) + smk] xf/k, 

k>Q 

where the operator = (c|,c_jt). In this form, it becomes 
apparent that the Ising model can be decomposed into a series 
of independent LZ transitions, as was first realized in [29]. 
Now, we can directly use the results of the previous section 
and write the supplementary Hamiltonian required for an adi- 
abatic driving across the QCP, 



^1 = -8 it) 



1 



sin^ 



k>0 



2 g^ + 1 — 2gcosk 



This expression is expected to be highly non-local in real 
space. It can be written as 



^1 = -g'{t) 



£ hm{gWP + -h^,2{gW, 
m=l 



[Nil] 



(6) 



The Hamiltonian ^^'"^ includes an interaction over a range 
m. Above, we have a factor of ^ for m — N /2 - which is 
the largest distance in the system with periodic boundary con- 
dition - because for even N there is only one spin over the 
distance N/2 while there are two different ones over smaller 
distances. Every is independent of g - all dependence 

on g is included in coefficients hm{g) - and reads 



The coefficients lim{g) are given by the Fourier transform, 
1 



hm{g) = -^/(fc)sin(m^), 



of the function 



sin^ 



4- g^ + 1 — 2gcosk 

In the limit of large we can approximate /?„ 
Jq f{k) sin{mk)dk/n, with the result 



(7) 



for 1^1 <1, 
for|g|>l. 
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Mapping back to spins, the supplementary Hamiltonians reads 



(8) 



Some comments are in order Firstly, since we have repre- 
sented the class of Hamiltonian in Eq. (4) - and in particular 
the Ising model (5) - as independent LZ crossings, the sup- 
plementary Hamiltonian ^ allows to adiabatically drive any 
eigenstate of the model under consideration [.>u]. This means 
that one can further tailor the Hamiltonian Jif{ for the purpose 
of driving exclusively a given subset of states, e.g., the ground 
states. Further considerations along that line of reasoning are 
beyond the scope of this letter [31]. 

Secondly, coefficients hm{g) can be neatly written as 
\h,nig)\ = e-('"±i)/^W/8, where ^(g) = |ln(|g|)|-' is the 
correlation length in the Ising model [32]. It follows that 
/^mfe) ~ ^(1) for distances m up to the correlation length and 
go to zero exponentially fast for a longer interaction range. At 
the critical point this means that Jifi is acting along the whole 
chain. 

Finally, Hamiltonians of the form J^^"'^ can be efficiently im- 
plemented in trapped-ion quantum simulators using strobo- 
scopic techniques [19, 33] already demonstrated in the labo- 
ratory [5]. 

Finite range interactions, filtering, and the KZM.- We next 
consider a linear quench of the coupling g{t) = gc — Vt, 
through the QCP at gc = 1, which brings the system from 
the paramagnetic to the ferromagnetic phase. The evolution 
induced by the Hamiltonian (6) is adiabatic in the instan- 
taneous eigenbasis of the Ising model (5) [m]. However the 
range of interaction, e.g. at the critical point, spans over the 
whole chain. As a result, from a practical point of view, one 
might be interested in assisting the crossing of the QCP with 
an approximation to Jfi that involves interactions of restricted 
range. In this section, we examine the simplest approxima- 
tion, namely, a direct truncation of Eq. (6) that hmits the range 
of interaction to M sites. 



M 

V S,„h,n 



(9) 



where s,„ will be a filter function. 

We start by examining the limit of a fast transition, which 
for now means D ^ I (this will be made more precise later). 
We consider the initial ground state in the paramagnetic phase 
and evolve it using: (i) only the supplementary Hamilto- 
nian = ^ (M), (ii) both the supplementary and the Ising 
Hamiltonians ^ ~ -j^o + ^ (M) . For both cases, we numeri- 
cally solve the time-dependent Bogoliubov - de Gennes equa- 
tions that describe the evolution of the system, as explained in 
[29]. Fig. 1 shows the probability of excitation in the A:-mode, 
Pk- 

Firstly, we have verified that pk does not depend on quench 
rate i) for I) ^ 1, and they coincide in both cases (i) and (ii) in 
that limit. That is, the presence of results only in phase 
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FIG. 1 . (Color online) Excitation probability as a function of the 
wave vector k, following a shortcut to adiabaticity in the \d quantum 
Ising model, (a) The crossing of the critical point is assisted by a 
truncated Hamiltonian .!^\ (M) {s^^ = 1) (b) and a modified truncation 
where the expansion coefficients are modulated by a raised-cosine 
Fourier filter. The range of the interaction increases from right to 
left with cutoff M = 4,8, 16,32,64 (N = 1600). Above, evolve the 
system from ^/ = 10 to = accros QCP at = 1 at quench rate 
D = 50 using = + (M) . The insets show the scaling of pi^ 
as function of kjkM ~ Mk. 



difference and does not affect how well the approximated 
supplementary-Hamiltonian is able to drive the transition. 
Secondly, the cutoff M in Eq. (9) implies approximating the 
function f(k) in Eq. (7) by its truncated Fourier series. Since 
f{k) is divergent and discontinuous at ^ = = 1 and A: = 0, 
M\ (M) suffers from so-called Gibbs phenomenon, this is, the 
problem of recovering point values of a nonperiodic or discon- 
tinuous function from its Fourier coefficients [34]. In Fig. 1(a) 
we present the results for the truncation Jifi (M) with s,n ~ 1 
(Dirichlet kernel). The Gibbs phenomenon is seen here in the 
appearance of side-lobes at large k. This can be prevented 
by using a Fourier space filter s,„ that modifies the expansion 
coefficients [34]. In Fig 1(b) we use the raised cosine filter 
■Sm = 5[1 + cos(m7r/M)]. It improves the convergence away 
from the discontinuity, making the decay of p^^ with k (almost) 
monotonic, and suppresses the side-lobes observed in its ab- 
sence at the expense of broadening p^. However, it remains 
impossible to recover f{k) close to its discontinuity, so modes 
with w are still excited. 

As an upshot, in the limit of fast quenches the effects of ap- 
proximating (M) by ^ (M) depends only on M, and to 
recover a fully adiabatic dynamics we need M = N/2. Based 
on the above considerations and the relation of /?„, to the corre- 
lation function in the Ising model, we can draw the conjecture 
that an approximation of the form in Eq. (9) induces an adi- 
abatic dynamics of the modes with k^ km ^ M^'. This is 
corroborated in the insets in Fig. 1 where we rescale k for dif- 
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FIG. 2. (Color online) Suppression of the total number of excitations 
«ex as a function of the quench rate u of the transverse field in a 
Id Ising chain following a quench through the QCP (g = 1). The dy- 
namics is assisted by a truncated auxiliary Hamiltonian Mi (M), with 
a cutoff M = 0, 1,2,4,8, 16,32,64 from top to bottom. The numer- 
ics for (m = 0) is the reference case where the KZM dictates a 
power-law scaling of iig^ for u ^ 1 . This scaling is recovered at slow 
quench rates in a passage through the QCP assisted by (M), while 
at faster rates there is an efficient suppression of excitations. As the 
range of the interactions is increased, the dynamics in all modes is 
driven through the instantaneous eigenbasis of Ji!o, and a complete 
suppression of excitations is achieved. N = 1600, we evolve from 
gi = 10 to = 0, and no filtering is applied (sm = 1). 

ferent values of the cutoff M, and the corresponding excitation 
probabilities collapse onto each other 
We consider as well the mean number of excitations, n^x = 
^ Jq Pkdk, as a function of quench rate v and cutoff M. The 
results are presented in Fig. 2. In a fast transition and for a 
given cutoff M, we are able to adiabatically drive modes with 
kM ^ ' . This means that the mean number of defects 
saturates at n^x ^ M-''. This limit can be seen at the right hand 
side of Fig. 2. 

Next, we focus on slower transitions which are induced 
by the Ising Hamiltonian and approximated supplementary 
Hamiltonian J^f — J^fo + -J/^i (M) . In the limiting case of 
M = the system is driven only by the Ising Hamiltonian and 
the non-adiabatic dynamics is correctly described by KZM 
[29, 35] - we do not consider here the super-adiabatic limit 
where the quench across QCP is adiabatic due to the gap re- 
sulting from the finite size of the system, which is reached for 
extremely slow quenches with v <^N^^ [12, 29]. 
For M = 0, there appears a characteristic value of momenta 
described by KZM: kKz ^ = v^^^ for the Ising 

model [30] (to be precise, we expect the power-law behavior 
for u <C 1, when the system goes out of equilibrium close to 
the QCP). The modes with k ^ kxz are expected to cross QCP 
adiabatically and = = v^^^. We recover this 

limit in the left hand side of Fig. 2, when M is small enough 
compared to i)-'. A crossover between the two quench rate 
limits occurs for intermediate values of M and v. The one 
which dominates for a given set of parameters depends on 



whether kKz is smaller or greater them km- 

Relation to the fidelity susceptibility.- Finally, it is interest- 
ing to draw a connection with the so-called fidelity suscepti- 
bility, which puts some constrains on ^ . Fidelity suscepti- 
bility, Xf(A), can be defined in the leading order expansion of 
the overlap of the ground-states calculated for slightly differ- 
ent values of external parameter A. For a finite system, in the 
limit 5 — > [37], we can Taylor expand the overlap in 5 [38]: 

\{X\X + S)\^^l-5hF{X). 

where Xfi^) for a non-degenerated ground state reads [30, 
39] 

^ mX)\d,M\n{X))\^ 

= — ^"^^oF — ' 

and {|«(A))} and E„ are instantaneous eigenstates and 
eigenenergies of Jf!o{X), and |0(A)) is the ground-state. 

In addition, the supplementary Hamiltonian J^i (2) which 
would be able to drive the evolution along the instantaneous 
ground state must satisfy 

(0(A)|^^KA))=a'(r)»l|i4ffi, (10) 

h„ — ho 

for n 0. Thus, we can verify that the mean variance of J^f 
in the instantaneous ground state is [23] 

Ajr2 = {0{X)\J^ff\OiX)) = |A'(r)|2xf (A), (11) 

Fidelity susceptibility, for translationally invariant system, is 
expected to typically scale [39-41] at the critical point as 
Xpi^c) ^ N^^''^, and away from the critical point as (A) ^ 
A^|A — Acl''*'^^. It is divergent in the vicinity of the QCP, as 
long as fidelity susceptibility is dominated by low lying exci- 
tations I |. 

In conclusion, for a broad family of many-body systems 
exhibiting a quantum phase transition, we have presented a 
method to assist the adiabatic crossing of the critical point at 
a finite-rate by supplementing the system with a finite-range 
time-dependent interaction. Our proposal is suited to access 
the ground state manifold in quantum simulators. The non- 
local terms of ,ifj''"'-type in the auxiliary Hamiltonian can be 
efficiently implemented using the stroboscopic techniques re- 
cently demonstrated in the laboratory [5, 19, 33]. We have fo- 
cused on the finite-rate adiabatic crossing of a quantum phase 
transition, where suppressing excitations is particularly chal- 
lenging due to the critical slowing down in the proximity of 
the critical point. Nonetheless, the method can be applied as 
well to the preparation of many -body states as an alternative to 
optimal control techniques [- :] or in combination with them. 
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